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Measurements of the evolution with redshift of the number density of massive galaxy clusters 
are used to constrain the energy density of massive neutrinos and so the sum of neutrino masses 
X] We consider a spatially-flat cosmological model with cosmological constant, cold dark matter, 
baryonic matter, and massive neutrinos. Accounting for the uncertainties in the measurements of 
the relevant cosmological parameters we obtain a limit of < 2.4 eV (95 % C.L.). 

PACS numbers: 14.60.Pq, 98.65.Dx, 98.80.Es 


Constraints on neutrino masses are of great interest 
for particle physics as well as for cosmology, and thus at¬ 
tract a lot of scientific attention (for recent reviews see 
Refs. 0, 0) S 0)- Current upper limits on the sum of 
neutrino masse^^ uii,, from cosmological structure for¬ 
mation data 0la, cosmic microwave background (CMB) 
fluctuation data HJaljOr combined CMB -|- large-scale 
structure data P- Tinlllli Elli , are of order an eV 0 
(for various limits see Table 1 of Ref. 0). The number of 
neutrino species can be constrained from Big Bang nu¬ 
cleosynthesis or by using CMB and large-scale structure 
data 00 

High energy physics experiments also constrain neu¬ 
trino masses, and have measured the number of light 
neutrino species with high precision 0. Direct searches 
for neutrino mass effects in beta decays yield limits in 
the region of several eV, but the sum over all neutrino 
masses is almost unconstrained by beta decay and other 
experiments, mainly due to the weak limit on the tau- 
neutrino mass. The measurement of neutrino oscillations 
on the other hand constrains the differences between the 
squared masses of the neutrino mass eigenstates Am^. 
With the justified assumption that neutrino masses are 
non-negative and for mass splitting Attiq rv 7 x 10“^ 
eV^ and Arn^tm > 1-3 ^ 10“^ eV^ 0, we obtain '^mi, 
> 0.04 eV if the solar mass splitting is between the high¬ 
est and second highest mass eigenstates (^ rrii, > 0.07 eV 
if the atmospheric mass splitting is between the two high¬ 
est states). Results from the LSND collaboration yield 
a larger lower limit on ^ m^, and must be considered if 
confirmed by the MiniBooNE experiment [l^, which is 
currently taking data. 

In this paper we use the dependence of galaxy 
cluster number density evolution on the massive neu- 
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trino energy density parameter Qi, to set a limit 
on We consider the standard spatially-flat 

ACDM Eriedmann-Lemaitre-Robertson-Walker cosmo¬ 
logical spacetime model with baryons, cold dark matter 
(CDM), massive neutrinos, and a non-zero cosmological 
constant A (for a recent review see Ref. 0 ) . To com¬ 
pute the cluster number density as a function of redshift z 
we use the Press-Schechter approach 00 as modified 
by Sheth and Tormen (ST) [2(]j.^ This approach makes 
use of the mass function N{M > Mq) of clusters (cluster 
number density as a function of cluster mass M greater 
than a fiducial mass Mn), which depends on cosmolog¬ 
ical model parameters [U, 0, 0, 0|. In particular, 
it is very sensitive to the matter density parameter UIm 
(= Ulb -I- Hcdm + where Ulb and Hcdm are the density 
parameters of baryons and CDM, respectively) and the 
value of (Tg (the r.m.s. amplitude of density fluctuations 
smoothed over a sphere of 8h“^Mpc radius, where h is 
Hubble constant in units of 100 km s“^Mpc“^) [0|. The 
observationally-viable ranges of these two parameters are 
related (for recent reviews see Refs. [13, |2a 1221 ), a 
current version of the relation between UIm and erg is 
given in Table 5 of Ref. 0| . The parameter erg is deter¬ 
mined by the matter fluctuation power spectrum which 
is sensitive to Uti, 00 . This is because the neutrinos are 
light particles and have a much larger free streaming path 
length^ than the CDM particles. Gravitational instabil¬ 
ity is therefore unable to confine the neutrinos on small 
and intermediate length scales, resulting in a suppres¬ 
sion of small- and intermediate-scale power. See Eig. 6 of 
Ref. 0 for cTg as a function of for models normalized 
to the Wilkinson Microwave Anisotropy Probe (WMAP) 


^ In the following we use PS to refer to the un-modified Press- 
Schechter approach. 

^ For ultrarelativistic particles the free streaming path length is 
equal to Hubble radius. After they became non-relativistic, par¬ 
ticle velocities are redshifted away adiabatically and so t he free 
streaming path length grows only slowly or decreases 
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data. 

Neutrinos are weakly interacting and this characterises 
how they affect cosmology. When rrii, > 10“^ eV neutri¬ 
nos are non-relativistic today 0 and thus behave like 
a hot component of dark matter. The presence of even 
a small fraction of massive neutrinos (hot dark matter) 
/y = riy/riM, of order 10-20 %, requires a smaller value 
for the cosmological constant in comparison to a pure 
ACDM model, while other cosmological parameters are 
largely unaffected @]. This is because a smaller IIa re¬ 
sults in a larger ^Im and hence a faster fluctuation growth 
rate, which compensates for the reduction of small- and 
intermediate-scale power caused by the neutrinos. Neu¬ 
trino free-streaming suppression of the linear^ growth of 
density perturbations on small and intermediate scales 
results in only a fraction of matter of order 1 — fi, be¬ 
ing involved in gravitational clustering |29| . Being an 
integral over the power spectrum, erg depends sensitively 
on This makes the cluster number density evolution 
with redshift very sensitive to the value of 
and so to the value of ^ rriu (since ^ rrii, = eV 

i) 

As with other cosmological tests, the cluster number 
density evolution test for neutrino masses requires fixing 
the range of some cosmological parameters. This may be 
viewed as a choice of priors; see Ref. 0 for a detailed dus- 
cussion of priors in the context of deriving neutrino mass 
limits from cosmological data. The parameter ranges 
we consider in this computation are picked as follows. 
Based on Hubble Space Telescope (HST) measurements 
of the Hubble constant 1^ . we use h = 0.71 ± 0.07 (one 
standard deviation limit).^ We assume adiabatic density 
perturbations with primordial power spectral index close 
to scale-invariant, n = 0.98 ± 0.02 (l-cr range) HHH 
Concerning the value of the matter density parameter 
Hm, there is evidence for S (0.2, 0.35) from different 
data such as Type la supernovae j WMAP , Sloan 
Digital Sky Survey (SDSS) and 2 degree Field Ga laxy 
Redshift Survey (2dFGRS) galaxy clustering WA . 
and galaxy cluster gas mass fraction evolution [2^. For 
a summary see Ref. who find 0.2 < ^ 0.35 

at two standard deviations; in our computation we use 
this as a one standard deviation range. We choose 
as ^ (0.77,1.11) as the two standard deviation range; 
for a discussion see Sec. 3.1 of Ref. [s^. For the baryon 
density parameter we use from Big Bang nucleosynthesis 


^ The effe ct of massive neutrino infall into CDM halos is studied in 
Ref. I28l . They found that three degenerate-mass neutrinos with 
52 ~ 2.7 eV alters the non-linear matter power spectrum by 

about 1 %. 

^ An analysis of all available measurements of the Hubble constant 
results in the more restrictive, but HST con sist ent, estimate h = 
0.68 di 0.07 (two standard deviation range) m. 


S (0.018, 0.022).® We work in a spatially-ffat model 
which requires 1 — Ha = + Hcdm -I- = Hm- The 

prior on Qb is less important than the priors on erg, Hm, 
h or n 1^. 

The paucity of galaxy cluster evolution data makes it 
inappropriate to use data analysis techniques based on 
fits. Reviewing such data-poor situations in the lit¬ 
erature one finds either modified fits assuming gaus- 
sian errors on the logarithm of the observed cluster num¬ 
ber density or likelihood approaches with Poisson errors 
based on the number of observed clusters. In our analy¬ 
sis we use a likelihood approach and we define our like¬ 
lihood function by £ = Hi Here ki is the 

number of observed clusters with mass greater that Mg, 
in the redshift bin (centered at redshift and of 
width Az) and /ii(M > Mg, Zi, Az) = Np,-ed{M > 

A/g, z) (iz/{Azai) is the predicted cluster number in this 
bin. Npred{M > Mo,z) is the predicted cluster number 
density. The normalization factor aj, which has dimen¬ 
sions of inverse volume, is the detection efficiency defined 
by NohsiM > Mo,Zi) = ah where Nohs{M > Mo,Zi) is 
the observed cluster number density in redshift bin. 

The predicted cluster number density Npred{M > 
Mo,z) depends on the cosmological model considered. 
An important characteristic of a cosmological model is 
the linear energy density perturbation power spectrum 
P{k,z). This is sensitive to the values of the cosmo¬ 
logical parameters h, Hm, and fib, as well as to the 
density parameter of each dark matter component, i.e., 
Hcdm and (although the requirement of flat spatial 
hyperspaces can be used to eliminate the dependence 
on Hcdm), and to additional parameters discussed below. 
This wavenumber-space two-point correlation function of 
density perturbations, P{k,z), at given redshift z, is de¬ 
termined by the power spectrum of initial perturbations 
Po{k), the energy density perturbation transfer function 
T{k,z), and the perturbation growth rate D{z). The 
functions T{k, z) and D{z) depend on model parameters 
and describe the evolution of density inhomogeneities 
[ifiji . To compute P{k,z) we assume a close to scale- 
invariant (Harrison-Peebles-Yu-Zeldovich) post-inflation 
energy density perturbation power spectrum Pg(fc) oc /c”, 
with n ~ 1, and we use a semi-analytical approxima¬ 
tion for the transfer function T(fc, z) in the AGDM model 
with three species of equal-mass neutrinos For the 


^ This is more consistent with estimates from WMAP data and 
from Big Bang nucleosynthesis using the mean of the primor¬ 
dial deuterium abundance measurements, but it is significantly 
larger than an estimate based on helium and lithium abundance 
measurements, see, e.g., Ref. m. 

® The effects of neutrino mass differences are irrelevant for our 
considerations, since if 52 mass eigenvalues 

are essentially degenerate 0, while if 52 <0.4 eV the mass 

differences do not much affect cluster number density evolution 
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growth rate D{z) we use Eqs. of Ref. |^, which 

are based on results from Refs. |42| . In summary, in our 
model the free parameters are n, h, ^Imj 
and (Tg (which fixes the normalization of the power spec¬ 
trum, as discussed below). 

The cluster mass function at redshift z is N{M > 
Mq,z) = J^^dM n{M,z), where n{M, z)dM is the co¬ 
moving number density of collapsed objects with mass 
lying in the interval {M,M + dM). In the PS approach 
the cluster mass function is determined by cr{R, z) the 
r.m.s. amplitude of density fluctuations smoothed over 
a sphere of radius R = (3M/AttpmV^^. where pM is 
the mean matter density [iq- The function n{M,z) is 
a universal function of the peak height Sc/cf{R), where 
Sc = 1.686. For gaussian fluctuations. 


n{M, z) oc 


a‘^{R,z) 


da{R, z) 
dM 


exp 


2a^{R, z) 


( 1 ) 


see, e.g., Eq. (1) of Ref. ^3- The evolution of the clus¬ 
ter mass function is determined by the z dependence of 
a{R,z). Now a‘^{R,z) is related to the power spectrum 
P{k,z) through 


a"(R,z) = ^ J Pik,z)\WikR)\‘^edk, (2) 
0 

where W(kR) is the Fourier transform of the top-hat win¬ 
dow function, W(x) = 3(sinx — x cosx)/x^. Numerical 
computation results for n{M, z) are not accura tely fit by 
the PS expression of Eq. (^1, see Refs. Sev¬ 

eral more accurate modifications of n{M, z) have been 
proposed, see Refs. US 0113. Here we use the ST 
modification 0, as defined in Eq. (5) of Ref. |24|. 


i{M, z) 


Sc 


a^{R, z) 


da{R, z) 


dM 


aSl 


<t‘^{R,z) 


-P 




\R,z) 


- 1/2 


exp 




2a‘^{R, z) 


(3) 


where the parameters a = 0.303 and p = 0.707 are fixed 
by fitting to the numerical results (for the PS case a = 1 
and p = 0) [001 • With this choice of parameter values 
the mass of collapsed objects in Eq. 0 must be defined 
using a fixed over-density contrast with respect to the 
background density pM |0 |2J, |0 [43 , and this requires 
accou nting fo r the mass conversion between Migoti and 
M 200 C [ 23 ,0."^ Such a conversion depends on cosmolog- 


^ The mass of a collapsed object is defined with respect to the 
Einstein-de Sitter critical density per = /{StvG), as the mass 
within a radius Rac^ inside of which the mean interior density is 
A times the critical density per- Assuming a p{R) density profile, 


ical parameters, see Fig. 1 of Ref. 0; we use an analyt¬ 
ical extrapolation of this figure to do the conversion for 
G (0.2,0.35). 

Our analysis is based on data from a compilation of 
massive clusters (with M > Mq = 8 x IO^'^M^Mq, where 
Mq here corresponds to the mass within a comoving ra¬ 
dius of 1.5/i“^Mpc and MqIs the solar mass) observed 
at redshifts up to z « 0.8, |0,0|, 


N{M > 8 X lO^^/i“^M0)/i-'^Mpc' 


X 10-^ 

1.7+fl X 10-®, 

1.4+i;i X 10-8, 
1.4+J f X 10-8, 


at 

at 

at 

at 


z = 0.0-0.1, 
z = 0.3-0.5, 
z = 0.5 — 0.65, 
z = 0.65-0.9. 


(4) 


Here the errors in N are from counting uncertainty and 
are automatically included in our Poisson likelihood anal¬ 
ysis. The bin width in z is chosen to reflect the uncer¬ 
tainty in the redshift measurement. 



p{R)R^dR = ApcrR%^c/^' -^ 200 c) which corresponds to 
an over-density A = 200 with respect to the critical d ensi t y, i s 
a common definition of the virial mass of the cluster l2li l24ll , 
while corresponds to the mass within a sphere of radius 

RiSOb inside which the mean density is 180 ti mes the background 
density pM {Risob = R54c for Qm = 0.3) Hi . For details of 
ma ss co nver sion assuming a Navarro-Frenk-White density profile 
m see m. 
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FIG. 1: The curves show the number density evolution of 
massive clusters (M > 8 x for models in which 

the sum of neutrino masses ^ m^=0, 2, 3.2 eV (from bottom 
to top). The other five parameters Qm, os, n, and h are 
set at the center of the scan interval. The crosses show the 
observational data of Eq. in four redshifts bins with I-ct 
P oisson error bars. 

We compute the likelihood from the data using Pois¬ 
son errors and the predicted number of clusters in each 
bin, and perform a maximum likelihood fit over a dis¬ 
cretized parameter space. We compare the observed clus¬ 
ter number density evolution of massive clusters with 
M > 8 X to model predictions for different 

values of ^ rrii, (see Figurenfor some examples) and for 
each value of ^ we marginalize the likelihood by inte¬ 
grating over the parameter space (0^, h, n, erg) with 
Gaussian weighting. Figure [3 shows the dependence of 
cluster number density evolution on ag. Here = 0 

and other four parameters {Q.m,^ b,h,n) are set at the 
center of the scan interval. 



FIG. 2: Variation of cluster number density evolution N{M > 
8 X 10^'^h~^M q, z) as a function of og with parameters kept 
at the same values as in Figure^except for erg which is varied 
and X) which is kept at zero. The observational data are 
the same as those shown in Figure 1. 

The likelihood as a function of ^ rrii, is shown in Figure 
13 Low values of ^ rrii, are favored but we find no pref¬ 
erence for a non-zero value. We obtain an upper limit of 
Y^rrii, < 2.37 eV (95 % C.L.), accounting for the uncer¬ 
tainties in all five cosmological parameters IIm, n, h, erg, 
and Ob. The cluster number density evolution depends 
signihcantly on erg (see Figure |3). The two-dimensional 
likelihood in (erg,^mjy) is shown in Figure 01 Lower 



FIG. 3: The marginalized (integrated) likelihood CjUmax as a 
function of X) on,,, accounting for the observational uncertain¬ 
ties in Om, Ob, erg, n, and h. The shaded region is excluded at 
95% confidence level. The likelihood function has been scaled 
for display purposes to a value of one at ^ m,, = 0. 


values of erg favor a higher neutrino mass. 

Our result of < 2.4 eV, which is based on the 

cluster number density evolution, is in good agreement 
with bounds on neutrino masses from CMB (and other) 
measurements 0,0 and corroborates that evidence. Our 
limit indicates that effects from the neutrinos on the evo¬ 
lution of galaxy clusters cannot be excluded and also in¬ 
dicates that these effects should be taken into account in 
the determination of cosmological parameters. Based on 
our result and on the particle physics limit of ^ rrii, > 
0.04 eV, we find Oj, is in the range of 0.1 % to 5 %. 
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